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We report on a combined theoretical and experimental study of the optical coupling between a 
microdisk resonator and a waveguide laying on different planes. While the lateral coupling between 
a planar resonator and a waveguide is characterized by a unique distance at which the resonant 
waveguide transmission vanishes because of destructive interference, the vertical coupling geometry 
exhibits an oscillatory behavior in the coupling amplitude as a function of the vertical gap. This 
effect manifests experimentally as oscillations in both the waveguide transmission and the mode 
quality factor. An analytical description based on coupled-mode theory and a two-port beam- 
splitter model of the waveguide-resonator coupling is developed, which compares successfully both 
to experimental data and numerical simulations. 
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In typical experiments with integrated circular-shaped 
resonators, known as Whispering-gallery mode (WGM) 
resonators, the cavity modes are probed by measuring 
the transmission spectrum of the bus waveguide [1] . The 
resonator-waveguide coupling (RWG) takes place via the 
overlap of the evanescent fields of the waveguide and res- 
onator optical modes. Therefore, the waveguide trans- 
mission spectrum shows a series of transmission dips at 
the cavity mode frequencies. In a planar geometry, i.e. 
when the waveguide and the resonator lay on the same 
plane, the mode coupling is known to grow exponentially 
when the gap x between these two is decreased. The 
parameters that are extracted from the waveguide trans- 
mission spectra are the transmission value at the bottom 
of each modal dip, T^jp, and its linewidth, F. The expo- 
nential variation of the coupling strength as a function of 
the gap X in this lateral coupling geometry results into 
clear trends in both Tdip and F. In particular, T^ip shows 
a single minimum at a critical gap x* at which the field 
propagating in the waveguide and the secondary emission 
by the resonator cancel each other on exact resonance due 
to a complete destructive interference For larger or 

smaller gaps, the transmission dips are shallower. On the 
other hand, following the exponentially growing RWC, 
the cavity mode linewidth F increases monotonically as 
the gap is reduced. 

So far, the optical coupling between a single-mode, rec- 
tilinear waveguide and a resonator has been usually de- 
scribed in terms of a spatially localized coupling at the 
minimum separation point. This provides an immediate 
access to the entire transmission spectrum, including key 
parameters such as the Tdip and the mode quality factor, 
Q, at the frequency ujq {Q — uq/T). The validity of this 
model in the lateral coupling geometry has proved itself 
in realizing in-plane all-integrated WGM resonator-strip 
waveguide systems with an optimized lateral gap that 



provides an efficient cross-coupling [3H7]- 

Recent experimental developments have shown that 
opting for a vertical coupling geometry for the RWC (the 
resonator and the buried bus waveguide lay on different 
planes) offers several advantages with respect to the lat- 
eral coupling geometry [5]: it does not require expensive 
lithographic techniques for gap definition and allows for 
an independent choice of materials and thicknesses for 
the optical components. While several experimental re- 
alizations of vertically coupled resonator/waveguide sys- 
tems have been reported [51-fTT]. no specific study of the 
influence of the vertical gap on the coupling efficiency has 
been performed yet. The issue of the vertical RWC has 
been discussed in Ref. 12J, where, the coupling efficiency 
optimization as a function of lateral [in-plane] alignment 
under fixed vertical gap conditions was addressed. 

In this Letter we show that, differently from the in- 
plane geometry, the critical coupling condition between 
a resonator and a bus waveguide in a vertical coupling 
geometry is fulfilled at several values of the gap distance 
X. This unexpected finding results from the spatially 
extended nature of the vertical RWC; the effective gap 
between the resonator and the waveguide remains almost 
constant for a sizeable distance Kf^ along the waveguide. 
Within this flat zone, the system can be modeled as a 
pair of coupled parallel waveguides. By using the cou- 
pled mode theory, we find that the interplay between 
the length A/z of the flat zone and the coupling length 
of parallel waveguides, Ac, results in an oscillating effec- 
tive coupling as a function of the gap x. This coupling 
is characterized by alternating regions of undercoupling 
and overcoupling separated by different critical coupling 
points. Our findings are in good agreement with exper- 
imental results and are further confirmed through full 
three dimensional (3D) numerical simulations. 

First, we start analyzing the geometry of an in-plane 
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FIG. 1. (Color online) 3D sketches of the geometry of (a) in- 
plane and (b) vertical coupled systems. In the in-plane case 
(a), the mode-to-mode separation (MM) corresponds to the 
distance L^or between the resonator and the waveguide. In 
the vertical case, MM is the distance L^er which is the hy- 
potenuse of the triangle formed by the two catheti, Lhor and 
X. As discussed in the text, the gap x is the distance between 
the two planes where the resonator and the waveguide lay. (c) 
The MM as a function of ^-coordinate along the waveguide 
for the in-plane (dotted line, Lhor) and the vertical (full line, 
Lver) geometries. In the vertical geometry, Lver exhibits an 
extended flat zone in the proximity of the minimum gap. The 
marked difference between the behavior of the first deriva- 
tives dLhor,ver/dz iu the two cases is shown in the inset. The 
resonator diameter is 50 /im. 



coupled system where, for simplicity, only the propaga- 
tion trajectories of the modes in an infinite strip waveg- 
uide and a WGM resonator, separated by a minimum 
gap of X are considered (Fig. [ija)). In this picture, 
the modes in the waveguide and the resonator propa- 
gate along the z-axis and a circle of a radius r, respec- 
tively. Then, the physical distance between the waveg- 
uide and the resonator (equivalent to a mode-to-mode 
separation) at a coordinate z (— r < z < r) is given by 
Lhor = X + r ~ Vr^ - z^- 

In a vertical coupling scheme, one deals with a 3D 
picture where the resonator and the waveguide modes 
lay on different planes, separated by a vertical distance 
X (Fig. [ijb)). Considering the simple case, when the 
resonator is vertically aligned to the z-axis at z — 
point, the mode-to-mode distance Ly^r reads as 
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At first sight, Lhor and L^er do not seem to be partic- 
ularly different. However, as Fig. [ijc) shows, Lyeriz) 



differs radically from Lhor{z), as the z runs from zero to 
±r. In fact, Lyer{z) is almost constant over a certain 
region around z = 0, which means that the variation in 
the waveguide-resonator separation here is insignificant. 
The extension A fz of the flat zone can be estimated from 
the inflection points of the second derivative d'^Lyer/dz^, 
as shown in the inset of Fig. [iJc); as expected, Ajz de- 
pends smoothly on x and is non-zero unless the vertical 
distance x vanishes. 

An important consequence of the existence of the flat 
zone is that the vertical RWC can be approximated to a 
coupling between two parallel waveguides within a good 
accuracy. This allows to employ the coupled-mode theory 
[l3l 114] to estimate the coupHng coefficients from the 
dimensions and material refractive indices of the devices 
used in actual experiments |15) . 

Figure [2] compares Af^ (for a 50/im diameter res- 
onator), and Ac (of two parallel waveguides) at a wave- 
length of A = 1.6 /im. The intersection of Afz{x) and 
Ac{x) has important physical implications. While Afz 
has a geometrical origin (resonator radius, vertical gap), 
the coupling length A^ depends on both the waveguide 
geometry and the refractive indices of the waveguides 
and the coupling medium. Depending of these, the expo- 
nential variation of Ac (a;) and the slowly varying Afz{x) 
across the crossing can lead to situations where: 

(A) Afz = Ac/2, which implies that half of the opti- 
cal power is transferred from the waveguide to the 
resonator according to cos^[7r2;/(2Ac)], 

(B) Afz^ Ac, which implies that all the optical power 
is transferred from the waveguide to the resonator, 

(C) Afz — 2Ac, which implies that all the optical power 
is coupled back into the waveguide by the end of the 
flat zone extension (Fig. [2]Ja),(b)). 

Based on these results, we proceed with modeling the 
amplitude of the transmitted light through the straight 
waveguide coupled to the WGM resonator as a function 
of frequency uj. According to the schematics shown in 
Fig. [2]jc), the resonator/waveguide coupling is modeled 
as a two-port beam splitter. For the output fields Eoi, 
Eo2 one has 
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The phase shift and the attenuation of the field after 
propagation along the rest of the cavity circumference 
is described by E^i = i;^! exp[j^(Lca^ - Afz)]. The 
refractive index of the resonator is considered to have 
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FIG. 2. (Color online) (a) The curves for the flat zone exten- 
sion (Ay 2, circles) and the coupling length (Ac, line) between 
parallel waveguides at A = show an intersection at 

a specific coupling gap. An important physical implication 
of the intersection is that critical coupling between the res- 
onator and the bus waveguide can manifest multiple times at 
different coupling gaps. A, B and C refer to the situations 
A/z = Ac/2, A/2 = Ac and A/z = 2Ac, respectively, (b) 
Coupled mode prediction cos^[7r2/(2Ac)] for the field inten- 
sity in the fiat zone z £ [0, A/z] in the three cases A, B and 
C. (c) Schematics of the vertical RWC reduced to an in-plane 
coupling model of a two-port beam-splitter. 



the form n = Jt-o + jf^j where no is the real effective 
refractive index, a is the overall intrinsic loss (material 
absorption, radiative and surface scattering) of the cavity 
with a circmxiference of Lcav = 27rr and A = l-Kcjix) is the 
wavelength. 

and 

— ^Lcav, the waveguide input and output fields are 



With the definition of two parameters 9 — 
■avj the waveguide input and output 
related as 



e-J"^ - cos 9 



and, thus, the waveguide transmission is T{u)) = 

\Eo2/Ei2\'^. 

After some algebraic transformation to Eq. |3j it is pos- 
sible to rewrite the waveguide transmission in the form 
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where Suj = ujq — w, a=cos9, b ~ — 
and 
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From Eq. |4j the typical Lorentzian lineshape of the 
cavity resonance, centered at luq and having a loss rate of 
r, can be recognized. In a weak coupling regime, when 
9 <^ 1, T tends to its lowest value Tq. The associated 
intrinsic (unloaded) Q-factor of the cavity is therefore 
Qo = ^^^^^ (for relatively weak intrinsic loss, ar <^ 1). 
Equations ^ and ([5| provide close expressions for the 
resonant waveguide transmission, T{uJo,x) and the Q- 
factor, Q{ujq,x), as function of the vertical gap x. These 
can be directly compared to experimental results. 

In our experiments, 50 /xm-diameter SiN^, planar 
WGM resonators were fabricated and coupled vertically 
to single-mode SiN^, waveguides as shown in Fig. [sj^a) 
(the only variable among different samples is the verti- 
cal gap thickness x). We refer to [8' for a more detailed 
discussion of fabrication and characterization techniques. 

From waveguide transmission experiments we have 
measured T{ujq, x) and Q{u!o, x) which are plotted as data 
points in Fig. [s] (b) and (c), respectively. Then, we ap- 
plied the developed theory to model the experimental re- 
sults of several first radial family modes at A ~ 1.59 /im. 
Figure [sj^b) shows that in fact, more than one minimum 
manifests in the waveguide transmission when the gap is 
monotonically reduced to zero. This result is in contrast 
with the in-plane case, where both the theory and ex- 
periments show a minimum at a unique critical coupling 
gap PHZj- However, this result is in agreement with the 
hypothesis of the "flat zone/coupling length" interplay. 
The computed Q{ujq, x), in its turn, shows multiple peaks 
with decreasing vertical distance (Fig. [sjjc)), while in the 
in-plane coupling scheme one observes a monotonic spoil- 
ing of the Q because of the exponentially increasing cou- 
pling loss. The computed values are compatible with the 
experimental data both for T and for Q-factor, though 
the experimental points are not numerous. It is interest- 
ing to note that the proposed model does not require an 
intermediate evaluation of the coupling constant and that 
by varying a one fits T{u;o,x) and Q(ljo,x) simultane- 
ously. Thus, by fitting the experimental data with Eqs.|4] 
and [5} an intrinsic loss of the resonator a ^ 2 cm^^ is es- 
timated, which gives a reasonable intrinsic Qq of 38,000. 

The agreement between the model and the data is even 
more evident when we compare the experimental data 
with the results of the in-plane model. We note here, 
that, within this model (Sj, T{u!q,x) and the Q{ujo,x) 
are indirectly related. In particular, using the exper- 
imental T{ujQ,x), one estimates the coupling constant 
K = Y^^J and, then, fits K with an exponential func- 
tion of the gap. This procedure requires an estimation 
of the critical coupling gap Xc for which if ~ 1, in or- 
der to switch between "=f" and "±" for the undercoupled 
and overcoupled regimes. Afterwards, the mode quality 
factor is calculated as Q{ujo,x) — Qo/{l + K{ojo,x)). In 
Figs. [3]jb) and (c) we report T(ljJo,x) and Q{ujo,x), ob- 
tained using the in-plane model, which however do not 
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troscopic analysis of the transmitted field, it suffers an 
intrinsic limitation of numerical methods, i.e. the errors 
can be determined by the finite resolution and time-steps 
of the performed simulations. The FDTD results are also 
reported in Fig. |3] (c) and show that at very small x the 
cavity modes have from two to three orders of magnitude 
larger Q's than those estimated from the in-plane model. 

In conclusion, we have shown that in a vertical cou- 
pling geometry between a WGM resonator and a bus 
waveguide, a series of critical coupling points appear at 
several values of the vertical gap. This phenomenon orig- 
inates from modal interference in the spatially extended 
coupling zone. By using a simple analytical model of the 
waveguide-resonator coupling, we have predicted an oscil- 
lating behavior between the undercoupling and overcou- 
pling regimes as a function of the separation gap, which 
manifests both in the contrast and the linewidth of the 
dips in the waveguide transmission spectrum. The model 
provides a good agreement with experimental data and is 
further confirmed by numerical simulations. Our findings 
reveal an important and unexpected difference between 
vertical and in-plane coupling in optical circuits and can 
impact the design of photonic devices that are presently 
of high interest for both fundamental science [TB] and for 
technological applications [TI] [T7] . 

The authors wish to thank M.R. Vanacharla and N. 
Prtljaga for stimulating discussions. This work has been 
supported by NAoMI-FUPAT. 
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FIG. 3. (Color online) (a) An experimental transmission spec- 
trum of a 50 fim WGM resonator, where some first radial fam- 
ily modes with different azimuthal mode numbers AI (labeled 
as are observed. Inset, optical image of the resonator 

and the buried waveguide, (b) The oscillatory behavior of 
the calculated (Eq. |4| and the measured peak transmission 
(points). The dashed line is the result of the classical in- 
plane coupling theory, (c) The calculated Q(a;)-factor (Eq.js] 
line), the measured Q-factor (full circles), the FDTD results 
(triangles) and the Q-factor obtained from the classical in- 
plane model (dashed line). The horizontal dotted line refers 
to half of the intrinsic Q-factor. 



provide a satisfactory fit to either experimental trends. 
Though the T(ujq,x) might agree with the experimental 
data, Q{uj),x) is markedly different. 

Finally, we have performed full 3D finite-difference 
time-domain (FDTD) simulations of the fabricated sys- 
tem. We note here, that due to the fact that the simula- 
tions run in time-domain, the eigenmode linewidths have 
been extracted employing a harmonic inversion method, 
which essentially consists in fitting the cavity ring-down 
signal with a trigonometric decay function. While this 
algorithm is much robust with respect to a simple spec- 
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